AD-A040  797  NAVAL  RESEARCH  LAB  WASHINGTON  D C 

APPLICATION  OF  SPEED  IMPROVEMENT  TECHNIQUES 
MAY  77  B H CANTRELL 
UNCLASSIFIED  NRL-MR-3S20 


END 

DATE 

FILMED 

' 7-77 


NRL  Memorandum  Report  3520 


f 

f . 

1 


?> 

o 

O 

/**\ 

< 

o 

<c 


Application  of  Speed  Improvement  Techniques 
in  Coordinate  Transformations  for 
the  Multiple  Platform  Sensor 
Integration  Problem 


B.H.  Canthf.ll 
Radar  Analysis  Stuff 

* Radar  Division 

: 


May  1977 


Q- 

O 


NAVAL  RESEARCH  LABORATORY 
Washington,  l).C. 


Approved  lor  public  release,  distribution  unlimited 


D DC 

rrsr^nn 
JUN  22  1977 
t^flbUU 


i 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (When  Dele  Entered) 

j’"~~  REPORT  DOCUMENTATION  PAGE 


M RE  RORT  NUMBER 


|2  GOVT  ACCESSION  NO. 


NRL\  Memorandum  Report  3520 


4.  I. rut  i*na  Suwiiti 

APPLICATION  OF  SPEED JMPROVEMENT  TECHNIQUES  IN 
, COORDINATE  TRANSFORMATIONS  FOR  THE  MULTIPLE 
; PLATFORM  SENSOR  INTEGRATION  PROBLEM” 


I? 


Ben  H.  Cantrell 


IS  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Naval  Research  Laboratory 
Washington,  D.C.  20375 

II  CONTROLLING  OFFICE  NAME  AND  ADORESS  12.  REPORT  DATE 

Department  of  the  Navy  May  1977 

Naval  Ship  Engineering  Center  ,,  NUMSER  of  paces 

Washington,  D.C.  20362  8 

"77  MONITORING  AGENCY  NAME  a ADDRESS (It  dllterenl  front  Controlling  Office)  15.  SECURITY  CLASS,  (ot  (him  report) 

UNCLASSIFIED 

15#  OCCL  ASSl  FIC  A TlON  DO W N&RA~D|  N G 
SCHEDULE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3 RECIPIENT'S  CATALOG  NUMBER 


5 TYPE  OF  REPORT  i PERIOD  COVEREO 

Interim  report  on  a continuing 
1 NRL  problem. 

[1  PERFORMING  ORG  REPORT  NUMBER 


• CONTRACT  OR  GRANt  NUM  B L Ri 

f)\  / ^ £ ■' 


10  PROGRAM  ELEMENT  PROJECT,  TASK 
AREA  * WORK  UNIT  NUMBERS 

NRL  Problem  R12-18 


M6  DISTRIBUTION  ST  ATEMEN  T <ol  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


jcc:ssiqh  tv 


17  DISTRIBUTION  ST  ATEMLNT  < ot  the  abatrecl  entered  In  Block  20,  If  different  from  Report)']  SmA:  t •tm 


**!!»  S«»ln  m 
Jlffjeflm. . q_| 


he  supplementary  notes 


19  KEY  WOROS  (Continue  on  reverse  aide  ll  nec  eatery  end  Identity  by  block  number) 

Command  and  control 
Radar  integration 
Computer  algorithms 

20  ABSTRACT  ( Continue  on  revet  me  aide  ll  neceaaery  and  Identity  by  block  number) 


j JM»3W»,1WUU31llir  CODES 
i — : fN — _ and.  or  SPECIAL- 


A means  of  improving  the  speed  of  computation  in  performing  coordinate  transformation 
used  in  multiple  platform  sensor  integration  problems  is  described.  The  method  is  based  on 
quaternions. 


DD  1473  EDITION  OF  1 NOV  65  IS  OBSOLETE  } 

S N " 10  2-0  14-  660  l 


SECURITY  CLASSIFICATION  OF  TmiS  PAGE  Pete  Mnteted 


/ 


CONTENTS 


INTRODUCTION 

EXAMPLE  

THE  METHOD  . 
SUMMARY  


iii 


APPLICATION  OF  SPEED  IMPROVEMENT  TECHNIQUES 
IN  COORDINATE  TRANSFORMATIONS 
FOR  THE  MULTIPLE  PLATFORM  SENSOR 
INTEGRATION  PROBLEM 


INTRODUCTION 

In  multiple  platform  sensor  integration  a large  number  of  coordinate 
transformations  are  required.  Therefore  it  is  necessary  to  improve  the 
speed  of  calculations  as  much  as  possible.  There  are  at  least  two 
fundamentally  different  starting  points  in  which  coordinate  rotations 
can  be  derived.  The  first  is  the  use  of  direction  cosine  matrices  which 
are  well  known,  published,  and  used  repeatedly  in  the  literature.  The 
second  means  is  to  derive  the  coordinate  rotations  by  using  quaternions. 
This  method  may  not  be  as  well  known.  Rather  than  develop  the  quater- 
nions an  example  is  first  shown  using  the  direction  cosine  and  then  the 
quaternion  approach.  Since  both  methods  yield  the  same  results  for  the 
same  initial  conditions  they  are  equivalent.  In  fact  identities  relate 
the  two  methods.  It  is  the  purpose  of  this  report  to  describe  to  those 
members  of  the  sensor  integration  community  (who  may  not  be  aware  of  it) 
how  these  identities  can  be  used  to  improve  the  computation  time  in 
coordinate  rotations. 


EXAMPLE 


A simple  rotation  of  the  X,Y  coordinates  into  the  X.',  H*  coordinates 
by  the  anqle  nt  is  shown  in  Pig.  1. 

Note:  Manuscript  submitted  May  10,  1977. 
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Fig.  1 — Demonstration  of  coordinate  rotation 
using  direction  cosines 


The  equations  are  well  known  and  are  as  follows: 


X = X cos  a + Y sin  a 


Y = -X  sin  a + Y cos  a 


The  same  transformation  in  quaternions  is 


2 (X  + Y tan  a/2) 


1 + i an  a/2 


(1) 

(2) 


(3) 


- Y - 2|X*Y  "/2)  tan  a/2  . 

1 + tan5  a/2 


(4) 


We  will  shown  (1)  and  (2)  are  equivalent  to  (3)  and  (4).  Expandiny 
(3)  and  (4)  one  obtains 
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(1  - tan  a/2) 
(1  + tan2  a/2) 


2 tan  a/2 
1 + tan2  a/2 


-2  tan  a/2  x + (1  - tan2  a/2)  y (6) 


1 + t an  a/2 


1 + tan  a/2 


By  using  the  trigometric  identities 


cos  a = 


1 - tan2  a/2 
1 + tan2  a/2 


sin  a 


2 tan  ' /2 
1 + tan2  “/2 


we  find  (1)  and  (2)  are  equivalent  to  (3)  and  (4). 

THE  METHOD 

In  the  example  used  the  quaternion  technique  would  be  faster  to 
compute  because  only  one  trigometric  function  need  be  computed.  In 
addition,  the  tangent  of  the  half  angle  is  also  an  easier  computation 
as  discussed  subsequently.  The  object  then  is  to  compute  the  tangent 
of  the  half  angle  and  through  the  use  of  trigometric  identities  to  re- 
arrange the  direction  cosine  equations  so  that  the  least  computation  is 
obtained. 

Consider  the  following  sequence  of  computations: 

b = tan  a/2  (9) 


D = 2/(1  + b ) 
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X 


D (X  + bY) 


X 


(11) 


Y"  = Y-bD (X+bY)  (12) 

for  example.  Equation  (11)  and  (12)  are  equivalent  to  either  (1)  and 
(2)  or  (3)  and  (4)  when  (9)  and  (10)  are  substituted  into  (11)  and  (12). 
The  calculation  requires  the  evaluation  of  one  trigometric  function, 

5 multiplies  or  divides,  and  4 additions  or  subtractions.  More  complex 
rotations  can  be  improved  in  the  computation  time  by  applying  similar 
principles.  The  principal  is  to  modify  the  sine  and  cosine  functions 
by  trigometric  identities  so  that  they  are  in  terms  of  the  tangent  of 
the  half  angle  and  then  rearrange  the  direction  cosine  equations  for 
the  least  computation. 

A few  comments  are  made  on  the  evaluation  of  the  tangent.  Since 
the  range  of  values  of  a is  from  -180°  to  180°  the  range  of  values  for 
the  half  angle  a/2  is  -90°  + 90°.  However  only  the  values  from  0°  to 
90°  need  be  computed  because  of  the  odd  symmetry.  This  range  of  values 
can  be  further  reduced  by  using  the  trigometric  identity. 

tan  a/2  = 1/ [tan  (90°  - a/2)]  (13) 

for  values  of  a/2  between  45°  and  90°.  Therefore,  the  only  values 
required  for  computation  are  those  of  a/2  between  0°  and  45°.  The 
computation  can  be  made  very  accurate  with  the  table  lookup  and  a 
linear  interpolation  between  points. 


SUMMARY 


Coordinate  rotations  can  be  derived  in  two  fundamentally  different 
ways.  The  quaternion  representation  has  a speed  advantage  to  the 
direction  cosine  representation.  It  was  the  purpose  of  this  report  to 
point  out  to  those  in  the  sensor  integration  community  this  fact.  The 
basic  principle  which  can  be  applied  is  to  find  the  tangent  of  the  half 
angle  and  rearrange  the  direction  cosine  equations  using  trigometric 
identities  to  obtain  the  higher  speed  forms.  The  technique  is  applica- 
ble to  either  digital  computers  or  high  speed  dedicated  digital  hardware 
where  extremely  high  speeds  are  desired. 


